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Abstract 

. We consider the homogenization of parabohc equations with large spatially- 

•/^ I dependent potentials modeled as Gaussian random fields. We derive the homog- 

enized equations in the limit of vanishing correlation length of the random po- 
. tential. We characterize the leading effect in the random fluctuations and show 

that their spatial moments converge in law to Gaussian random variables. Both 
1^ ' results hold for sufficiently small times and in sufficiently large spatial dimensions 

■ d > m, where m is the order of the spatial pseudo-differential operator in the 

parabolic equation. In dimension d < m, the solution to the parabolic equation is 
shown to converge to the (non-deterministic) solution of a stochastic equation in 
the companion paper [2] . The results are then extended to cover the case of long 
range random potentials, which generate larger, but still asymptotically Gaussian, 
random fluctuations. 
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^ . 1 Introduction 

Let m > and P{D) the pseudo-differential operator with symbol p{^) = |.^|™. We 
consider the following evolution equation in dimension d > m: 

(^- + P{D)--q{^)y,{t,x) = 0, xGM^ t>0, 

Here, uo e L^(M'^) and q{x) is a mean zero stationary Gaussian process defined on a 
probability space JF, P). We assume that q{x) has bounded and integrable correlation 
function R{x) = ¥,{q[y)q{x + y)}, where E is the mathematical expectation associated 
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with P, and bounded, continuous in the vicinity of 0, and integrable power spectrum 
(27r)'^^(0 = /jjd e"'^ ''-R(x)(ix in the sense that J^d\B{o i) ^(Ol^l'^^^^^ < The size of 
the potential is constructed so that the hmiting solution as e — > is different from the 
unperturbed solution obtained by setting q = 0. The appropriate size of the potential 
is given by 

^.^plln.l^ d = ra, 

The potential is bounded P-a.s. on bounded domains but is unbounded P-a.s. on 
M*^. By using a method based on the Duhamel expansion, we nonetheless obtain that 
for a sufficiently small time T > 0, the above equation admits a weak solution ■) G 
L'^{il X W'-) uniformly in time t E (0, T) and < e < Eq. 

Moreover, as e — > 0, the solution converges strongly in L'^{Q x M.'^) uniformly 
in t e (0,T) to its limit u{t) solution of the following homogenized evolution equation 

(^ + P(Z})-p)M(t,x) = 0, xeM^ t>o, 

m(0, x) = Uo{x), X G 

where the effective (non-negative) potential is given by 



CdR{0) d = m, 



pd 




' ' -(24 d > m. 

Here, q is the volume of the unit sphere S'^~^. We denote by Qj^ the propagator for the 
above equation, which to uo{x) associates Q^uo{x) = u(t,x) solution of ([3]). 

We assume that the non-negative (by Bochner's theorem) power spectrum R{C,) is 
bounded by /(|^|), where /(r) is a positive, bounded, radially symmetric, and integrable 
function in the sense that r'^~^""^f{r)dr < oo. Then we have the following result. 

Theorem 1 There exists a time T = T{f) > such that for allt G (0,T), there exists 
a solution u^it) G L^(r2 x M*^) uniformly in < s < Eq. Moreover, let us assume that 
R{^) is of class C'(M'^) for some < 7 < 2 and let u{t,x) be the unique solution in 
L'^iW^) to ([3]). Then, we have the convergence results 

\\{Ue -Us){t)\\L2^n^Rd) < e2\\uo\\L^Rd), 

(5) 

||(u^ - u)(t)||i2(iRd) < e^^^||no||L2( 



where a < b means a < Cb for some C > 0, a A b = min(a,6), where Ue(t,-) is a 
deterministic function in L^(M'^) uniformly in time, and where we have defined 



llnel ^ d = m, 

^d-m m < < 2m, 

5"| lne| = 2m, ^ ' 

e'^ d> 2m. 



The Fourier transform ^^(t,^^) of the deterministic function Us{t,x) is determined ex- 
plicitly in below. 
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Note that the effective potential — p is non-positive. The theorem is vahd for times T 
such that 4Tpj < 1, where p/ is defined in lemma [2^ below by replacing R{C) by /(|^|) 
in the definition of p in (jlj). 

The error term — u is dominated by deterministic components when e"'^^ ^ 
and by random fluctuations when e'''^^ <^ In both situations, the random 

fluctuations may be estimated as follows. We show that 

Ul,e{t, X) = {Ue - E{Ue}) jt, x), (7) 

e 2 

converges weakly in space and in distribution to a Gaussian random variable. More 
precisely, we have 

Theorem 2 Let M be a test function such that its Fourier transform M & {W^) fl 
L'^iW^). Then we find that for all t e {0,T) 

Ke(t,-),M)-^ / Mt{x)adW,, Mtix)= [ gPM{x)glMx)ds, (8) 

JK.'i Jo 

where convergence holds in the sense of distributions, dW^ is the standard multiparam- 
eter Wiener measure on M*^ and a is the standard deviation defined by 

:= (27r)'^^(0) = / E{q{0)q{x)}dx. (9) 

This shows that the fluctuations of the solution are asymptotically given by a Gaussian 
random variable, which is consistent with the central limit theorem. 

We observe a sharp transition in the behavior of at = m. For d < m, the 
following holds. The size of the potential that generates an order 0(1) perturbation is 
now given by (see the last inequality in lemma 12.21) 

d 

= £2. 

Using the same methods as for the case d > m, we may obtain that ^^(t) is uniformly 
bounded and thus converges weakly in L^(fi x M"') for sufficiently small times to a 
function u{t). The problem is addressed in [2], where it is shown that u{t) is the 
solution to the stochastic partial differential equation in Stratonovich form 

^ + P{D)u + uoa^ = 0, (10) 
ot dx 

with u{0,x) = Uo{x) and ^ d-parameter spatial white noise "density". The above 
equation admits a unique solution that belongs to L^{Q x M'') locally uniformly in time. 
Stochastic equations have also been analyzed in the case where d >m (i.e., d >2 when 
P{D) = —A), see [9], [12]. However, our results show that such solutions cannot be 
obtained as a limit in L'^{Q x M^) of solutions corresponding to vanishing correlation 
length so that their physical justification is more delicate. In the case d = 1 and m = 2 
with q{x) a bounded potential, we refer the reader to [13] for more details on the above 
stochastic equation. 
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The above theorems [T] and [2] assume short range correlations for the random poten- 
tial. Mathematically, this is modeled by an integrable correlation function, or equiv- 
alently a bounded value for R{0). Longer range correlations may be modeled by un- 
bounded power spectra in the vicinity of the origin, for instance by assuming that 
R{0 ~ ^iO^iO^ where S{C,) is bounded in the vicinity of the origin and h{^) is a 
homogeneous function of degree — n for some n > 0. Provided that > m + n so that 
p defined in (|4]) is still bounded, the results of theorems [1] and [2] may be extended to 
the case of long range fluctuations. We refer the reader to theorem [3] in section 13.31 
below for the details. The salient features of the latter result is that the convergence 
properties stated in theorem [T] still hold with /3 replaced by /3 — n and that the random 
fluctuations are now asymptotically Gaussian processes of amplitude of order e 2 . 
Moreover, they may conveniently be written as stochastic integrals with respect to some 
multiparameter fractional Brownian motion in place of the Wiener measure appearing 
in (IH]). 

Let us also mention that all the result stated here extend to the Schrodinger equation, 
where ^ is replaced by i-^ in We then verify that —p in ([3]) is replaced by p so that 
the homogenized equation is given by 

[i-+p{D)+py{t,x) = 0. 

The main effect of the randomness is therefore a phase shift of the quantum waves as 
they propagate through the random medium. Because the semigroup associated to the 
free evolution of quantum waves does not damp high frequencies as efficiently as for the 
parabolic equation ([1]), some additional regularity assumptions on the initial condition 
are necessary to obtain the limiting behaviors described in theorems [T] and [2J We do 
not consider the case of the Schrodinger equation further here. 

The rest of the paper is structured as follows. Section [2] recasts ([1]) as an infinite 
Duhamel series of integrals in the Fourier domain. The cross-correlations of the terms 
appearing in the series are analyzed by calculating moments of Gaussian variables and 
estimating the contributions of graphs similar to those introduced in [3 [H]. These 
estimates allow us to construct a solution to ([T]) in L^(f2 x M.'^) uniformly in time for 
sufficiently small times t G (0, T). The maximal time T of validity of the theory depends 
on the power spectrum R{C,)- The estimates on the graphs are then used in section [3] to 
characterize the limit and the leading random fluctuations of the solution Ue{t,x). The 
extension of the results to long range correlations is presented in section 13. 3[ 

The analysis of ([T]) and of similar operators has been performed for smaller potentials 
than those given in ([2]) in e.g. [H [B] when converges strongly to the solution of the 
unperturbed equation (with g = 0). The results presented in this paper may thus be seen 
as generalizations to the case of sufficiently strong potentials so that the unperturbed 
solution is no longer a good approximation of n^. The analysis presented below is 
based on simple estimates for the Feynman diagrams corresponding to Gaussian random 
potentials and does not extend to other potentials such as Poisson point potentials, let 
alone potentials satisfying some mild mixing conditions. Extension to other potentials 
would require more sophisticated estimates of the graphs than those presented here or 
a different functional setting than the x R'') setting considered here. For related 

estimates on the graphs appearing in Duhamel expansion, we refer the reader to e.g. 

mum- 
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2 Duhamel expansion and existence theory 

Since q{x) is a stationary mean zero Gaussian random field, it admits the following 
spectral representation 

^ ^ e^^-^Q{dO, (11) 



q[x) 



where Q{dC,) is the complex spectral process such that 



for all / and g in L^(M'^; R{^)d^) with the power spectrum and correlation function of 
q respectively defined by 

0<{2itYR{^)= [ e-'^''R{x)dx, Rix) = E{q{y)q{x + y)}. (12) 
In the sequel, we write Q{dO = q{Od^ so that E{g(^)g(C)} = R{0^i^ + C) and 

2.1 Duhamel expansion 

Let us introduce (jeiO — ^'^""^(^O? ^^e Fourier transform of e~"q'(|). We may now 
recast the parabolic equation ([T]) as 

{^^+C)Ue = qe*Us, (13) 

with u^(0,C,) = Uo{^), where 

qe*%{t,0= I Ue{t,i-OQe{dO= [ Ue{t, ^ - C)UC)dC ■ 

Here and below, we use the notation ^"^ = l^l"". After integration in time, the above 
equation becomes 



^,(t,0 = e-*«"no(0+ fe-'^'" [ q^{^ - ^,)u^{t - s,C,)d^,ds. (14) 

Jo Jr'' 

This allows us to write the formal Duhamel expansion 

^.(^,0 = 5^^n,.(t,0, (15) 

riGN 

<e{t, Co) = / n / e-«"^'=e-(*-S^^o Yl g,(a - ^k+i)MCn)dsde (16) 
Here, we have introduced the following notation: 

n— 1 n 

s = (so, . . . , Sn-i), tfc(s) = t - So - ... - Sk-1, to{s) = t, ds = Y\ dsk, d$, = Y\ d^k- 

k=0 k=l 



We now show that for sufficiently small times, the expansion (fT5|) converges (uni- 
formly for all e sufficiently small) in the L^(f2 x W^) sense. Moreover, the norm of 
is bounded by the L'^{W^) norm of Mq, which gives us an a priori estimate for the 
solution. The convergence results are based on the analysis of the following moments 



f/;'"^(t, C) = E{ne,n(t, i)Ue,m{t. C)}, (17) 

which, thanks to (ITB]) . are given by 

-1 ptk{s)m-l ^t,i-r) 



r '111 rk(s) '"-^ /-tix-T) 

n—l m—1 



fc=0 1=0 

Let us introduce the notation s„(s) = i:„(s) = t — X]fc=o -^fc ^^"^ Tm(''") = tm{T) = 
t - Ya=o ^i- ^Iso define ^n+k+i = Cm-k and for < k < m. Since 

is real-valued, we find that 

/n+m+l n+m 
H e-^^«"E{ n g.(a-a+i)}^o(en)4(en+i)cisrf|, 
fc=0 k=0,k^n 

where the domain of integration in the s and ^ variables is inherited from the previous 
expression. Note that no integration is performed in the variables s„(s) and s„_|_i(t). 
The integral may be recast as 

/n+m+l n+m n n+m+l 

H e-^'=«'E{ H g,(e,-efc+i)}Mo(en)^o(a+i)5(t-$^Sfc)5(t- Yl ^k)dsdt 

k=0 k=0,k=/=n k=0 k=n+l 

where the integrals in all the Sk variables for Q<k<n + m + l are performed over 
(0, oo). The 5 functions ensure that the integration is equivalent to the one presented 
above. The latter form is used in the proof of lemma 12.11 below. 

We need to introduce additional notation. The moments of u^^n are defined as 

f/;(t,e) = EK,„(t,o}. (18) 

We also introduce the following covariance function 



K"'™(t, C) = cov(n,,„(t, 0, UeAt. 0) = u^'^^it, e, C) - u^it, Omt^ C). (19) 

These terms allow us to analyze the convergence properties of the solution Ue(t, ^). Let 
M(^) be a smooth (integrable and square integrable is sufficient) test function on M*^. 
We introduce the two random variables 



hit) = / \ueit,0\d^ (20) 
Xe{t) = [ n,(t,0^(Oc?e (21) 
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2.2 Summation over graphs 



We now need to estimate moments of the Gaussian process q^. The expectation in 
jjn,m. Yamshe^ unless there is n G N such that n + m = 2n is even. The expectation 
of a product of Gaussian variables has an explicit structure written as a sum over all 
possible products of pairs of indices of the form — ^k+i- The moments are thus given 
as a sum of products of the expectation of pairs of terms qei^^k — ^fe+i), where the sum 
runs over all possible pairings. We define the pair (^fc, ^;), 1 < /c < /, as the contribution 
in the product given by 

mMk^i - ik)Uii~i - ii)} = - ik-i)mk - 6-1 + - ez-O- 

We have used here the fact that -R(— = -^(0- 

The number of pairings in a product of n + m = 2n terms (i.e., the number of 
allocations of the set {1, . . . , 2n\ into n unordered pairs) is equal to 

(2.-1). (^^ 

There is consequently a very large number of terms appearing in ?7"''^(t, ^o, ^n+m+i)- 
In each instance of the pairings, we have n terms k and n terms / = l{k). Note that 
l{k) >k + l. We denote by simple pairs the pairs such that l{k) = /c + 1, which thus 
involve a delta function of the form 5{^k+i — ■Cfe-i)- 



Figure 1: Graph with n = 3 and m = 1 corresponding to the pairs (^1,^3) and {^2,^5) 
and the delta functions 6{^i — .^0 + {3 ~ (,2) and 6{C,2 — G + ■Cs ~ ^4)- 

The collection of pairs {C,k,^i{k)) for n values of k and n values of l{k) constitutes a 
graph 6 constructed as follows; see Fig {1] and [5]. The upper part of the graph 
with n bullets represents u^^^ while the lower part with m bullets represents Ue,m- 
The two squares on the left of the graph represent the variables ^0 and C,n+m+i in 
[/"'"^(t, ^0, ^n+m+i) while the squares on the right represent 'Uo(^n) and uo{^n+i)- The 
dotted pairing lines represent the pairs of the graph g. Here, (J5 denotes the collection 
of all possible \&\ = 2^-i{fi-iy. gi'aphs that can be constructed for a given n. 

We denote by Aq = Ao{g) the collection of the n values of k and by Bq = -Bo(g) the 
collection of the n values of l{k). We then find that 

n+m 

k=0,k^n ge© k&Ao(g) 
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This provides us with an exphcit expression for U^'^{t, C,o, Cn+m+i) as a summation over 
all possible graphs generated by moments of Gaussian random variables. We need to 
introduce several classes of graphs. 

We say that the graph has a crossing if there is a. k < n such that l{k) > n + 2. 
We denote by 0c C the set of graphs with at least one crossing and by &nc = 0\0c 
the non-crossing graphs. We observe that V^"'"^(t, ^n+m+i) is the sum over the 
crossing graphs and that f/"(t, ^o)U"^{t, ^n+m+i) is the sum over the non-crossing graphs 

in[/-'-(t,eo,en+™,+i). 

The unique graph with only simple pairs is called the simple graph and we 
define (£>ns = ^\Qs- We denote by &cs the crossing simple graphs with only simple 
pairs except for exactly one crossing. The complement of 0cs in the crossing graphs is 
denoted by &cns = ^A^cs- 

As we shall see, only the simple graph Qs contributes an 0(1) term in the limit 
£ — > and only the graphs in (3cs contribute to the leading order 0{€2^'^~'^"^) in the 
fluctuations of u^. 

The graphs are defined similarly in the calculation of U^(t, ^o) in (|T8|) for n = 2n and 
m = 0, except that crossing graphs have no meaning in such a context. A summation 
over k G ^0(0) of all the arguments — ■Cfc-i + ^i{k) — 6(fc)-i of the 6 functions shows 
that the last delta function may be replaced without modifying the integral in U^{t,C,o) 
by 5(^0 -U- 

This allows us to summarize the above calculations as follows: 

/n+m+l 
k=0 ge© (22) 

n £"-'"^(£(6 - 6-i))5(6 - 6-1 + - ^m-i)dsdt 

fceAo(g) 

Similarly, 

/n 
k=o ge© (23) 

n ^'-'"i?(^(a - 6-i))5(a - 6-1 + - ^m-i)dsdt 

fe6Ao(g) 

2.3 Analysis of crossing graphs 

We now analyze the infiuence of the crossing graphs on /e(t) and X^(t) defined in (1201) 
and (l^T]) . respectively, for sufficiently small times. We obtain from (1191) and (1221) that 

/n+m+l 
gtv^c k=o (24) 
n £'^-'°^(£(efc - 6-i))5(6 - 6-1 + ^i(k) - ^m^i) ds 

involves the summation over the crossing graphs <3c. Let us consider a graph g G C5c 
with M crossing pairs, M > 1. Crossing pairs are defined hj k < n and l{k) > n + 2. 
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Denote by {C,q„-,,^i{q^)), 1 < m < M the crossing pairs and define Q = maxmiqm}- By 
summing the arguments inside the delta functions for all k < n, we observe that the 
last of these delta functions may be replaced by 

M 
m=l 

Similarly, by summing over all pairs with k > n + 2, we obtain that the last of these 
delta functions may be replaced by 

M 
m=l 

The product of the latter two delta functions is then equivalent to 

A/-1 



m=l 



The analysis of the contributions of the crossing graphs is slightly different for the energy 
in fl20l) and for the spatial moments in fl2Tl) . We start with the energy. 

Analysis of the crossing terms in /e(t). We evaluate the expression for | V^"'"*(t, ,^0; ■Co) 
in fl2^ at C,n+m+i = and integrate in the .^o variable over M'^. Let us define A' = 
Ao\{Q}. For each k G A' U {0}, we perform the change of variables —* We then 
define 

k^A'u{o} 

^ keA'u{0}. 

Note that = Cn+i since ^n+m+i = Co- This allows us to obtain that 



&=\t:iz:'.zi (25) 



» » n+m 

■^^^ geaJc k=i 

n e-''^m-ea-Mj~Ck-i + ^m-^kk)-i) (26) 

feGA'(0) 

M 

e-"^mo -e^n + Yl - ^q„_i)5(en,+i - in)dsdt 

m=l 

Here ci^ also includes the integration in the variable .^o- The estimates for V^'^ here 
and in subsequent sections rely on integrating selected time variables. All estimates are 
performed as the following lemma indicates. 

Lemma 2.1 Let t > given and consider an integral of the form 

n— 1 ftkXs) 

= n / ( n ^'^(^^)) n (2?) 

k=0 k=0 k=0 
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where < /^(s) < 1 for < k < n and assume that < h At. Then 

4-2 < (/^At)/„_i. (28) 

Moreover, let s be a permutation of the indices < k < n — 1. Define I^-i o-s In-i with 
fk replaced by fs{k)- Then = /„_i. 

Using the above result with the permutation leaving all indices fixed except s{n — l) = 
K and s{K) = n — 1 for some < K < n — 2 allows us to estimate by integrating 
in the Kth variable. 

Proof. The derivation of ( l28i) is immediate. We also calculate 

n— 1 n n 

k=0 k=0 k=0 

71—1 n n 

( n Uk){s,ik)))5{* - Yi n ^^fc 

fc=0 fc=0 fc=0 

71—1 n n 

fc=0 fc=0 fc=0 

□ 

Note that e"'^"'-^''^^"-*" and e~'^"+^''^^^^= ^™ are bounded by 1. We now estimate the in- 
tegrals in the variables sq, Sn+m+i, and Sk for A; G A' in (l26l) . Note that n + 1 cannot 
belong to A' and that n does not belong to A' either since either n = Q (last crossing) 
or n G -Bo is a receiving end of the pairing line k — > l{k). Each integral is bounded by: 



' \ — 

m 



ds < — At. (29) 

— Cm ^ ' 



The remaining exponential terms e ^^^^k) are bounded by 1. Using lemma [2711 this 
allows us to obtain that 



At) 

fceA'(0) ^k 



2 



771=1 



Here, ds corresponds to the integration in the remaining time variables Sk for k ^ 
A' U {n + m + 1}. There are 2?2 — 1 — (n + 1) = n — 2 such variables. Note the square 
on the last line, which comes from integrating in both variables sq and Sn+m+i- 

The delta functions allow us to integrate in the variables C,i(k) for k G A'{g) and 
the initial condition uo{^n) in the variable Thanks to lemma [2l2] below, the power 
spectra allow us to integrate in the remaining variables in A'U {0}. The integrals in the 
variables in A' are all bounded by defined in lemma fl?2\ whereas the integral in ^0 
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results in a bound equal to e^p/, where is defined in (|6]). As a consequence, we have 
the bound 



/ |K"'™(t,eo,eo)Meo < E ( / d^pY'\\uo\?Pfe^ = E ( / d^p"} 



Using Stirling's formula, we find that |0c| < 2"-"(n-i)! bounded by (^)". It remains 
to evaluate the integrals in time. We verify that 

11/ dsQ- ■ ■dsn-\ = tk{s) = t - So - . . . - Sk-i- (30) 

Let p = p{g) be the number of Sk for /c < n in s and g = q{g) be the number of Sk for 
/c > n + 1 in s, with p + q = n — 1. Using fl30l) . we thus find that 



rfs = -- = < <r-^fi 



p\ q\ {n — p J \ 2e J \2e 

using Stirling's formula. This shows that 

E ( /rfs) < ^(4p;T)^ (31) 



uniformly for t G (0,T). We thus need to choose T sufficiently small so that 4p/T < 1. 
Then, for r such that 4p/T < < 1, we find that 

J lK"''"(^,eo,eo)Meo < Cr"+'"£^||nof , (32) 

for some positive constant C. It remains to sum over n and m to obtain that 

|E{/,(t)} - / E{u,{t,OVd^\ < -.e^Wuof. (33) 

We shall analyze the non-crossing terms generating |E{'Ue(t, ^)}p shortly. Before doing 
so, we analyze the infiuence of the crossing terms on X^. We can verify that the error 
term in fl33|) is optimal, for instance by looking at the contribution of the graph with 
n = m = 1. 

Analysis of the crossing terms in X^. It turns out that the contribution of the 
crossing terms is smaller for the moment than it is for the energy I^. More precisely, 
we show that the smallest contribution to the variance of X^ is of order e'^-^a f^j, graphs 
in and of order 5'^-2a+/3 for the other crossing graphs. 

We come back to (124|) and this time perform the change of variables ^ — for 
k & A' only. We re-define 

^>^-\ k keA\ ^^^^ 
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and find that 

/n+m+l 
W e-^'=(«^)™no(en)^o(en+i) 

gt^Dc fc=0 

n e-'-m - eiU)5{^j - + - ^^^^ 

e''-^^R{e{iQ - e^_i))5(e„+™+i - a+i + in- io)dsdi. 

Note that neither n nor n + m + l belong to A'{q). For each /c G ^'(fl), we integrate in 
Sfc and obtain using ( 129|) that 

ge^c"^ fc^A'(g) 

n (7^ ^) ^^^^ - - + - ^m-i) ^^^^ 

keA'(s) 

By assumption on -R(i^), we know the existence of a constant i?oo such that 

e''-'"R{e{^Q-eQ^,))<e'"'-R^. (37) 

This is where the factor e'^-^o ^^iggg^ \Yg need however to ensure that the integral in 
is well-defined. We have two possible scenarios: either Q = n 01 n E Bq. When Q = n, 
the integration in is an integration in for which we use Uo{^n)- When n G Bq, 
we thus have n = l{ko) for some ko and we replace the delta function involving by a 
delta function involving given equivalently by 

M 

- ^Q-i + ^o-^n + Y,U- U-i). (38) 

m=l 

In either scenario, we can integrate in the variable C,q without using the term R{e{C,Q — 
^Q_i)). We use the inequality 

\uo{^n)Uo{^n+l)\ < ^{iM^n)]'^ + |Mo(^n " ^0 + ^n+m+l) , (39) 

to obtain the bound 

\K''''^{t,^o,Um+i)\ < e'-'"Roo J2 ( / ds)p'}-'\\uof. (40) 

The bound is uniform in ,^0 and ^n+m+i- Using fl3Tl) and fl32l) . we obtain 

lK"''"(t,eo,en+m+i)| < £'^"'"r"+™||nof . (41) 
After summation in n, m G N, we thus find that 

E{(X, - E{XJ)2} < ^£^-2"||Moiri|M||?. (42) 
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Similarly, by setting ^n+m+i = ^o, we find that 

e{ /" \us\\t,oviOdd- [ iEK(t,o}ivm <^^'-'ii^oirii¥^iii, (43) 

for any test function (f G //^(M*^). This local energy estimate is to be compared with 
the global estimate obtained in (!33|) . 

Analysis of the leading crossing terms in X^. The preceding estimate on may 
be refined as only the crossing graphs in C5cs have contributions of order e'^"^". We 
return to the bound (136!) and obtain that 

\V,^'"'{t,^o,Um+i)\<e''-^^RooY, [ n e-^''(«'')"|no(e„)^o(en+i)| 

n ^"'"(t^ ^ ^)r(^>^ - ^^i-Mj - ^i-i + - kk)-i) ^^^^ 

fc6A'(g) 

+m+l 

The n + 3 variables in time left are Sq, Sn+i, sq, Si(Q), and the n — 1 variables s;(A'(g)). 

Let G (5c- Let us assume that for some k such that {^k, 6(fc)) is not a crossing pair, 
we have l{k) — 1 > k, i.e., g G (S„cs- The non-crossing pairs are not affected by the 
possible change of a delta function involving ^„ to a delta function involving ^q. We 
may then integrate in the variable si(^k) and obtain the bound for the integral 

— + — ^o) 



thanks to lemma 12.21 below. The summation over all graphs in <Sncs of any quantity 
derived from VJ^''^{t,C,o,^n+m+i) is therefore smaller than the corresponding sum 
over all graphs in 0^. We thus see that any non-crossing pair has to be of the form 
l{k) — 1 = k, i.e., a simple pair, in order for the graph to correspond to a contribution 
of order 

Let us consider the graphs composed of crossings and simple pairs. We may delete 
the simple pairs from the graph since they contribute integrals of order 0(1) thanks to 
lemma 12.21 below and assume that the graph is composed of crossings only, thus with 
n = m and Q = n after deletion of the simple pairs. Let us consider k < n with 
l{k) > n + 1 so that the delta function 

- Ck~l + kk) - ^i(fe)-l) 

is present in the integral defining V^'"^. We find for the same reason as above that 
the contribution of the corresponding graph is of order e'^~^°'e^ by integration in the 
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variable si^k)- As a consequence, the only graph composed exclusively of crossing pairs 
that generates a contribution of order e'^-^o jg ^^le graph with n = m = 1. This concludes 
our proof that the contribution of order e'^-^o? yn,m jg given by the nm graphs in C5cs 
when both n and m are even numbers (otherwise, is empty). All other graphs in 
0c provide a contribution of order smaller than what we obtained in (HTj) . In other 
words, let us define 

/n+m+l 
gtCc. fc=o (45) 

fceAo(0) 
We have found that 

|K"'"^(t,eo,en+,n+i) - K:r(t,eo, W+OI < e'^-^^+V+^Hnoir. (46) 
2.4 Analysis of non-crossing graphs 

We now apply the estimates obtained in the preceding section to the analysis of the 
moments U^{t) defined in (fT8|) and given more explicitly in (!23|) . Our objective is to 
show that only the simple graph g contributes a term of order 0(1) in (!23|) whereas all 
other graphs in C5„s contribute (summable in n) terms of order 0{e^). Note that n = 2n, 
for otherwise, U^{t) = 0. We recall that the simple graph is defined by l{k) = k + 1. 
We thus define the simple graph contribution as 

?7"(t,eo) = w;(t,eo)«o(eo) 



Kit.^o) = /ne-^^^«™n\'^-2"i2(£(6fe+i-6fe))5(6{fc+i)-6fc)c?srf^, ^^^^ 

k=0 k=0 



and 



UeAt.^o) = 5^f/".(i,eo) := W.(t,eo)%(eo). (4 

nSN 

For all k G Aq, we perform the change of variables ^ ^ and (re-)define as before 

^k k^Ao 



10- 



This gives 



/n 
n 



e 



060 fc=o ^50) 
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Assuming that l{k) — 1 > k for one of the pairings, we obtain as in the analysis leading 
to (H6|l the following bound for the corresponding graph: 



This shows that 
so that 



1 



\E{u,}{t,0 - UeAt,0\ < ^^^l^o(OI, (52) 

at least for sufficiently small times t G (0, T) such that 4pfT < 1. It remains to analyze 
the limit of Us^s{t,C,) to obtain the limiting behavior of and I^^ip- This analysis is 
carried out in the next section. Another application of lemma [2^2] shows that U^^siii is 
square integrable and that its L^(M'^) norm is bounded by \\uo\\- In other words, we have 
constructed a weak solution Usit) ^ L'^{fl x M'^) to ( fT3l) since the series ( fT5l) converges 
uniformly in x W^) for sufficiently small times t G (0,T) such that 4p/T < 1. 

Collecting the results obtained in (133!) and (!52|) . we have shown that 

IKUe - f/e,s)(t)||L2(nxMd) < £ 2 || % || L2(Kd) , 

where ^/g^^ is the deterministic term given in (1481) . The analysis of f4 and that of is 
postponed to section [HI after we state and prove lemma 12.21 which allows us to analyze 
the contributions of the different graphs. 

Lemma 2.2 Let us assume that R is bounded by a smooth radially symmetric, decreas- 
ing function f{r). We also assume that f{r) < Tfr^" for some 0<n<d — xnin 
dimension d > m and n = when d <vn.. Then we obtain the following estimates. 
For d > m, we have 

/I - f°° I 

^R(^, - y)d^, < := c, ^/(l^l)l^l'^-irflel V r;, 

uniformly in y ^ W^, where q = |S''^~^| and a\/ b = max{a, b). Moreover, 



For d = m, we define pf = q/(0) and have 



lne\ 1 = 1 
1 = 2. 



For d <m, we have 
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Proof. Once R is bounded above by a decreasing, radially symmetric, function /(r), 
the above integrals are maximal when y = z = thanks to lemma 12.31 below since 
|^|~"^ and (e"^!^!""^ A t) are radially symmetric and decreasing. The first bound is then 
obvious and defines pj. The second bound is obvious in dimension d > 2m since I'Cfcl"^'" 
is integrable. 

All the bounds in the lemma are thus obtained from a bound for 



We obtain that the above integral restricted to r G (l,oo) is bounded by a constant 
times e^pf for d >vn. and by a constant times for d < m. It thus remains to bound 
the integral on r G (0, 1), which is equal to 



Replacing /(r) by r/-r~", we find that the first integral is bounded by a constant times 
e'^~" and the second integral by a constant times e'^~^ V e'"^ when d — n — Im. ^ and 
e^'"! lne| when d = 2m — n. It remains to divide through by when / = 2 to obtain 
the desired results. □ 

Lemma 2.3 Let f , g, and h be non negative, hounded, integrable, and radially sym- 
metric functions on M'^ that are decreasing as a function of radius. Then the integral 



which is well defined, is maximal at ( = t = 0. 

Proof. In a first step, we rotate ( to align it with r. The first claim is that the 
integral cannot increase while doing so. Then we send ( and r to 0. The second claim 
is that the integral again does not increase. 

We assume that the functions /, g, and h are smooth and obtain the result in the 
general case by density. We choose a system of coordinates so that r = |r|ei, where 
(ei, . . . , ed) is an orthonormal basis of M'^, and ( = \(\9 with 9 = {cos 9, sin 9, 0, . . . , 0). 
Without loss of generality, we may assume that 9 G (0,7r). Then I(^^r may be recast as 
le and we find that 



where we denote = h{0 with the same convention for / and g and define 






(53) 





It is sufficient to show that dgJg < 0. We find 



doJf 



e 



L 



Sd-i 



9^-Vfmi^-09mi'-r)d^, 
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with 6-^ = (— sin 9, cos 9,0, ... ,0). We decompose the sphere a.s ip = {ip ■ 9, ip) and find, 
for some positive weight w{fi) that 

We now observe that 

as 11^1(6' ■ 11)9 + ip) — t\ < \\S,\{9 ■ il)9 — ip) — t\ by construction. Indeed, we find that 
\\^\{9-^lj9±^)-T\'^-\i\'^-\T\^ + 2\T\\e,\9-ij9-T = ±2|r||^|^-r = ±2|r| |^|^^ ■ r whereas 
9^ ■ T = — sin^|r| < by construction. This shows that \^\{9 ■ il)9 + ip) is closer to r 
than 1^1(6' ■ ip9 — ip) is, and since g{r) is decreasing, that deJg < 0. This concludes the 
proof of the first claim. 

If /? = or r = 0, we set 6 = below. Otherwise, we may assume without loss of 
generality that r = —hC, for some h > 1. We still define C, = \C\9. We now define the 
integral la = lacx^ < a < 1, and compute 

dala= [ -C-^f{^-aOg{^ + bOh{Od^= [ -(.VfiOgi^+ib-aKM^ + aQd^. 

Define /(^, Q = + (b — a)()h{^ + a(). Then because / is radially symmetric, we have 

daL= mi\m\'~'d\^l Mlel) = -/'(lel) / 9.^^;m^,C)d^|J. 
Jo Js'i-^ 

We recast 

mm) = -rm) f (o ■ - i{-mx))d^ < o, 

since | + 7C| > | — + 7C| by construction for all 7 > and thus for 7 = and 
J = b — a. This shows that dala < and concludes the proof of the second claim. □ 

3 Homogenized limit and Gaussian fluctuations 

In this section, we conclude the proof of theorems [1] and [2j 
3.1 Homogenization theory for 

We come back to the analysis of Us,s(t,^) defined in fl47l) . Since only the simple graph 
is retained in the definition of mean field solution Ue^s{t,0^ ^^e equation it satisfies 
may be obtained from that for by simply assuming the mean field approximation 
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KlqsqeUs} ~ ^{QeQey^iue} since the Duhamel expansions then agree. As a consequence, 
we find that U^^s is the solution to the following integral equation 

ft f ij~ — s I* 

t V 

= e-*«™no(0+ / r e-ri^-s^)^~^?s,,d-2a f _ ^^^^^^^^^ - v,Od^ids,dv 

Jo Jo ^ ^, J 

Jo Jo J 

(54) 

The last integral results from the change of variables eC,i — > and sie — >■ si. It 
remains to analyze the convergence properties of the solution to the latter integral 
equation. Note that ^ acts as a parameter in that equation. Let us decompose 

AeU{t, = Ps f e-^'^-Uit - V, Odv + E,U{t, 0, (55) 
Jo 

with = Jjgd ^^^^m'^^^ d^i when d > m and = CdR{e^) when d = m. Then we have 

Lemma 3.1 Let ,^ G M'^ and f{r) as in lemma [27B . Then the operator E,, defined above 
in (l55l) is bounded in the Banach space of continuous functions on (0,T). Moreover, we 
have 

||^.||£(C(0,T)) <£^-" (56) 

Proof. We start with the case d > m so that and e^^'^a _ Note that n in lemma 
12.21 is defined such that c/ > m — n as well. With = — in (l55i) . we find that 



BeUUt^i) = / e-«"^ / / e-^'^'^R{i,-ei)d^rds,UU-t-v.Odv. 
Jo Jo J 

The remainder E'g is then given by 

E.UUt.i) = f r I e^^'"''ie'''^"''^-l)e-^^'^m,-eOd^ids,U,,sit-v,Odv 
Jo, Jo J 



t POO 



e-« ''e-«i ''R{^i - e^d^idsiU.At - i)dv. 

The continuity of E^Ui;^s{t,0 time is clear when [/^^^(t,,^) is continuous in time. 
Without loss of generality, we assume that Us^si'^O is bounded by 1 in the uniform 
norm. We decompose the integral in the Si variable in the first term of the definition of 
E^ into two integrals on < Si < ^ and ^ < si < Because e~^"^(e'^"'^"'^^ — 1) < 1, 
the second integral is estimated as 

r / e-«"-(e^"«"^^ - l)e-«"^^i?(ei - eOd^ids.dv 
Jo J^ J 

fj ^e~^^^m-emidv< j ^{^^M)m-emi<e'~''Pf, 
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< 



thanks to lemma [2?2l The above bound is uniform in ^. The last integral defining on 
the interval Si > ^ is treated in the exact same way and also provides a contribution 
of order 0{e^~'^). 

The final contribution involves the integration over the interval < Si < Using 

m mm 

g-5 f (^gE i ^1 — 1) < e'^^'^sie 2~ on that interval, it is bounded by 
Jo Jo Jr'' ^ 

n mem r 

< ^(1 - e-V) / / e-^"-i?(6 - ei)di,ds,, 

Jo JR'i 

by switching the variables < s < < Using lemma 12.3^ we may replace 

— ^0 by -R(^i) in the above expression. This shows that 



JM'* Jo 



Jo 

We observe that 

2 



so that 



Jo SI 



poo , 



The integral over (1, 00) is bounded by e'^p/. Using the assumption that /(r) < r 
we obtain that the integral over (0, 1) is bounded by a constant times 



r m p\ 



r'^-^-"dr + / ^rf-i-n-2m^r < r^-^ V 1, 

/I ~ ' 

Jr^m 

when d — n — 2m 7^ and | lnr| when d = n + 2m. Since r is bounded by a constant 
times e~'^, this shows that J3 is bounded by e^^-"^-^ when (i — n — 2m 7^ and e*^! lne| 
when d = n + 2m. This concludes the proof when d > xn — n. 

We now consider the proof when d = m with n = 0. Then, _ -j^^. The 

leading term is given by [4 s, which solves the integral equation: 



1 



lln^l Jo Jo 

(57) 

Here we have defined 

B,Uit,0=Pe[ e-^''''Uit-s,Ods, p, = c,R{eO, 
Jo 
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and Ef, is the remainder. As in the case c? > m, a contribution to | lne|i?£ comes from 
Jo Jo J 

We again decompose the integral in si into < Si < 2^ and < Si < We have 

t ^ 

I r I e-«"'^(e^"«"^i - l)e-«"^ii?(ei - eOd^idsidv 
Jo J^ J 

/2 /p"^ \ 2 ^ 

according to lemma 12.21 Also, 

£ I e-«"^(e^"^"- - l)e-^"-^(ei - ^O^^i^^i^^ < e-{^ V l) 

according to the calculations performed above on J3, which is uniformly bounded, and 
thus provides a | Inel"^ contribution to E^. 
We are thus left with the analysis of 

U{t, ^ j - ei)di, - p,) U{t - V, Odv, 

as an operator in £(C(0,T)) for fixed. Define Re{C,i) = R{C,i — ^O- The integral in ,^1 
may be recast as 

1 -e-T" 



R,{re)di2{e))dr. 

We observe that the integral on (1, 00) is bounded by ||-R||i. Assuming that R is of class 
C°'^(M'^) for 7 > 0, we write ^^(^1) = ^^(0) + (^^(Ci) - ^^(0)). The second contribution 
generates a term proportional to r'^ in the integral and thus is bounded independent of 
e. It remains to estimate 



1 



CdRM / dr = CdRsiO) / dr. 

Jo ^ Jo 

The latter integral restricted to (0,1) is bounded. On r > 1, e~''"/r is uniformly 
integrable so that 

CdReiO) / ^— rfr = CdR{eO\lne\ + 0(1). 

Jo r 

This shows that E^ is of order = as an operator on C(0, T) and concludes the 

proof of the lemma. □ 

Note that may be written as 



Jo 
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where (ps{s,C,) is uniformly bounded in s, ^, and e by a constant v^oo- The equation 

{I-A,)U{t,0 = S{t,0, 

admits a unique (by Gronwall's lemma) solution given by the Duhamel expansion and 
bounded by 

\u{t,o\<\\s\W^-. 

As in the proof of lemma ISTTl let us define -B^ = — E^. We verify that 
the solution to 

(/-i?,)it, = e-*«"n,(0, 

is given by 

H.(t,0 = e-*(«"-''^^^"^o(0- (58) 
The solution may thus grow exponentially in time for low frequencies. The error 
Veit,0 = (t/,,,(t,e) -il,(t,0) is a solution to 

iI-A,)V, = E,ii,{t,0, 

so that over bounded intervals in time (with a constant growing exponentially with time 
but independent of ^), we find that 

\Veit,0\<e'. (59) 
Up to an order 0{e'^\uo{^)\), we have thus obtained that K{ue{t,^)} is given by 

which in the physical domain gives rise to a possibly non-local equation. It remains to 
analyze the limit of the above term, and thus the error PeiO ~ Pi which depends on the 
regularity of R{i). For R{£) of class C2(M'^), we find that 



|g-t{r-p.{0) _ e-t(r-p)i 



The reason for the second order accuracy is that -R(— = -^(0 and V-R(O) = so 
that first-order terms in the Taylor expansion vanish. For R{C) of class C"'{W^) with 
< 7 < 2, we obtain by interpolation that 

|e-t(r-P.«)) _ e-i(r-P)| < e^*e-«"V^tC. 

When m > 7, the above term is bounded by 0{e"') uniformly in ^ and uniformly in time 
on bounded intervals. When m < 7, the above term is bounded by 0{e'^) uniformly in 
^ and uniformly in time on bounded intervals. This concludes the proof of theorem [TJ 
In terms of the propagators defined in (H7|) . we may recast the above result as 

|W,(t,0-W(t,0| <^"^^ ZY(t,0 = e-«" (60) 
where the bound is uniform in time for t G (0, T) and uniform in ^ G W^. 
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3.2 Fluctuation theory for 

We now address the proof of theorem [2j The first term in the decomposition of n„ ^ 
defined in fll6l) is its mean E{n„£}, which was analyzed in the preceding section. The 
second contribution corresponds to the graphs 0cs in the analysis of the correlation 
function and is constructed as follows. Let n = 2p + 1, p G N. We introduce the 
corrector ^ given by 

n p 



/II, tj ^ 

fc=0 g=0 r=l 

r ^ 1 



r=q+l 



In other words, all the random terms are averaged as simple pairs except for one term. 
There are p + 1 such graphs. We define 



(^'0 = E<,^(^'^)- (62) 



n>l 



We verify that 

is equal to the sum in V^'"^{t, ^o, ^n+m+i) only over the graphs in (3cs- Indeed, the above 
correlation involves all the graphs composed of simple pairs with a single crossing. 
Now let us define the variable 

Y, = {ue-ut-nUe},M)- (63) 

Summing over n,m G N the inequality in ( 146|) as we did to obtain (l42l) . we have 
demonstrated that 

^{Ye} <e''-^'"^^\\uof\\M\\l (64) 

for sufficiently small times. The leading term in the random fluctuations of Us is thus 
given by u'^. It remains to analyze the convergence properties of 

ZM = ^iK,M). (65) 

e 2 

We thus come back to the analysis of and observe that for n = 2p + 1, 

p ^ q 



/II, f H 

fc=0 g=0 r=l 

r ^ 1 

Ui0-^n)[ n ^'"'"%(6r-en,))5(e2r.-l-a)Jwo(en)c?srfC- 



r=q+l 
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Using the propagator defined in fHTj) . we verify tliat 

, 2g+l 



« — n 7 — n ^ — 1 



g=0 ^ ^-=0 r=l 

ft 



V / w,'^(t-t2,+i,eo)g.(eo-uwr'''(^2,+i,e«)%(e«)rft2,+ic?e« 



q=0 

Upon summing over n, we obtain 

utit, = [ m - s, - ^iHis, ^i)uo{^i)dsd^i. (66) 

We can use the error on the propagator obtained in (160|) to show that the leading order 
of is not modified by replacing by U. In other words, replacing Ue by U modifies 
Ze in (1651) by a term of order 0{e^^'^^'^^) in L^{il x W^), which thus goes to in law. 

Note that u'^{t,C,) is a mean zero Gaussian random variable. It is therefore sufficient 
to analyze the convergence of its variance in order to capture the convergent random 
variable for each t and C,- The same is true for the random variable Z^. Up to a 
lower-order term, which does not modify the final convergence, we thus have that 

{ut,M) = j j\^,{t-S,me{il)Uu,{s,i-il)dsdidi,. 

We have defined Uf(t,^) = W(t,^)/(^) for a function f{^). As a consequence, we find 
that, still up a vanishing contribution, 

° X w,^(s,e-a)Wno(r,c-Ci)c?[5rcciai]- 

Here and below, we use the notation d[xi . . . Xn] = dxi . . . dxn- By the dominated 
Lebesgue convergence theorem, we obtain in the limit 

E{|Z|2}:=i?(0) j I j j\,{t-s,i)Uu,{s,i-ii)diids'\di. 

Here, Z is defined as a mean zero Gaussian random variable with the above variance. 
Let us define Gtf{x), the solution at time t of ([3]) with f{x) as initial conditions, which 
is also the inverse Fourier transform of Uj{t,^). We then recognize in J J^l^Mit — 
s,^)Uuo{s,^ — ^i)d^ids the Fourier transform of J^t{x) defined in so that by an 
application of the Plancherel identity, we find that 

E{Z^} = {27iyR{0) [ ([ gj:^,M{x)gPuo{x)dsydx={27iyR{0)[ M^t{x)dx. (67) 

This shows that Z(t) is indeed the Gaussian random variable written on the right hand 
side in (|H]) by an application of the Ito isometry formula. This concludes the proof of 
theorem [21 
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3.3 Long range correlations and correctors 

Let us now assume that 

m = Hosio, < h{xo = ixrno, m 

where h{^) is thus a positive function homogeneous of degree — n and S{^) is bounded 
on -8(0, 1). We assume that -R(0 is still bounded on M'^\-B(0, 1). We also assume that 
m + n < d and that p in (j4]) is still defined. We denote by (f{x) the inverse Fourier 
transform of h{$,). Then we have the following result. 

Theorem 3 Let us assume that h{C,) = \^\~" for n > and m + n < d. We also impose 
the following regularity on uq: 

[ |{to(e + r)\'^h{^)d^ < C, for all t E W^. (69) 

Jb{o,i) 

Then theoremU\ holds with (3 replaced by j3 — n. 
Let us define the random corrector 

uiAt,x) = -J^iue-E{u,}){t,x). (70) 

e 2 

Then its spatial moments {ui^£{t, x), M{x)) converge in law to centered Gaussian random 
variables Af{0,T,M(t)) with variance given by 

SM(t) = (27r)'^^(0) / MtixMx - y)Mt{y)dxdy. (71) 

Proof. The proof of theorem [1] relies on three estimates: those of lemma 12.21 and 
lemma 13.11 and the uniform bound in (1371) for R. Lemmas 12.21 and 13.11 were written to 
account for power spectra bounded by in the vicinity of the origin. It thus remains 
to replace ( 1371) by 

when I^Q — ^Q_i\ < 1 while we still use ( 137|) otherwise. We have defined 6*00 as the 
supremum of S'(^) in -6(0, 1). It now remains to show that the integration with respect 
to in (l36l) is still well-defined. Note that either Q = n or — Cq-i may be written 
as — C for some C G M'' thanks to (l38l) . Upon using (l39l) . we thus observe that in all 
cases, the integration with respect to $,q in (136|) is well-defined and bounded uniformly 
provided that (1701) is satisfied uniformly in r. Using the Holder inequality, we verify 
that (1701) holds e.g. when Uo{- — r) G L'^{B{0, 1)) uniformly in r for g > This 
concludes the proof of the first part of the theorem. 
Let us now define 

ZM = -^{u^,,M)=e'^Z,{t). 

e 2 

We verify as for the derivation of E{Z^} that 

¥.{Zl} = 1 1' l'u^^^{t-s,Ol^M{t-r,OS{e^imimi-Ci) 
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The dominated Lebesgue convergence theorem yields in the hmit e ^ 



where Ait is defined in (jH]). An apphcation of the inverse Fourier transform yields ( 17T1) . 
□ 

Note that fITT]) generalizes fl67|) . where (p{x) = S{x), to functions Mt{x) E L^(R'^) with 
inner product 

{f,g)v= f{x)g{y)v{x - y)dxdy. (72) 

For = 1^1"", we find that ip{x) = c„|x|"~^, with c„ = r(^)/(2"7rir(f )) a nor- 
malizing constant. Following e.g. [7l[TU], we may then define a stochastic integral with 
fractional Brownian 

Z = [ Mt{x)dB^{x), (73) 



where is fractional Brownian motion such that 

E{B^{x)B^{y)} = l{\xr + \yr -\^-yn, 2H=1 + ^. 

We then verify that E{Z^} = E^/ so that the random variable Z is indeed given by the 
above formula fl73l) . When n = 0, we retrieve the value for the Hurst parameter H = ^ 
so that B^ = W. The above isotropic fractional Brownian motion is often replaced in 
the analysis of stochastic equations by a more Cartesian friendly fractional Brownian 
motion defined by 

d 

^h{x) = l[H,{2H,-l)\x,\^''^-\ 

1=1 

The above is then defined as the Fourier transform of 

d d 

i=l 1=1 ^ 

The results of theorem [1] and [3] may also be extended to this framework by modifying 
the proofs in lemmas 12.21 and 13. 1[ We then obtain that (1731) holds for a multiparameter 
anisotropic fractional Brownian motion S^, H = {Hi, . . . , Hd), with covariance 



E{B^{x)B^{y)} = 1 n d^^l'''' + l^^l'""' - 1^' - y'^l"""')- 



1=1 



Note that homogenization theory is valid as soon as c? > m + n. As in the case n = 0, we 
expect that when d < xn + n (rather than d < m), the limit for will be the solutions 
in L^(f2 X M'^) to a stochastic differential equation of the form (ITOj) with white noise 
replaced by some fractional Brownian motion; see also [8]. 
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The stochastic representation in fl73l) is not necessary since (t) fully characterizes 
the random variable Z . However, the representation emphasizes the following conclu- 
sion. Let and Z^ be the limiting random variables corresponding to two moments 
with weights M\[x) and M^ix) and a given Hurst parameter H . When H = |, we deduce 

directly from (175]) that K{Z^Z2} = when Mi{x)M2{x) = 0, i.e., when the supports 
of the moments are disjoint. This is not the case when H ^ ^ a.s fractional Brownian 
motion does not have independent increments. Rather, we find that E{Zf^Z|^} is given 
by {-Mt,!, Mt,2)ip, where the inner product is defined in fl72|) and Ait,k is defined in dH]) 
with M replaced by Mk, k = 1,2. Similar results were obtained in the context of the 
one-dimensional homogenization with long-range diffusion coefficients [3]. 
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